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One of the most exciting subjects in solid state physics is a single layer of graphite which exhibits 
a variety of unconventional novel properties. The key feature of its electronic structure are linear 
dispersive bands which cross in a single point at the Fermi energy. This so-called Dirac cone is 
closely related to the surface states of the recently discovered topological insulators. The ternary 
compounds, such as LiAuSe and KHgSb with a honeycomb structure of their Au-Se and Hg-Sb 
layers feature band inversion very similar to HgTe which is a strong precondition for existence of the 
topological surface states. In contrast to graphene with two Dirac cones at K and K' points, these 
materials exhibit the surface states formed by only a single Dirac cone at the V point together with 
the small direct band gap opened by a strong spin-orbit coupling (SOC) in the bulk. These materials 
are centro-symmetric, therefore, it is possible to determine the parity of their wave functions, and 
hence, their topological character. Surprisingly, the compound KHgSb with the strong SOC is 
topologically trivial, whereas LiAuSe is found to be a topological non-trivial insulator. 
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The search for new materials with inverted band struc- 
ture provides the basis for the Quantum Spin Hall effect 
(QSH)[l|-[l3|. This new exciting field started with the 
prediction and experimental observation of the QSH in 
quantum wells in two-dimensional topological insulator 
of the binary semiconductor HgTe [l|, Q . In a series of 
single crystals, such as Bi2Se3, three-dimensional topo- 
logical insulting behavior was observed in topological sur- 
face states appearing as Dirac cones 7 9]. Later on, the 
manifold of Heusler semiconductors with 18 valence elec- 
trons and a similar band inversion was proposed [lol [llj . 
Since this proposed class of materials is extremely rich, 
it provides much wider flexibility in design by tuning the 
band gap size and the spin-orbit coupling (SOC) mag- 
nitude. In addition, the multi-functionality allows the 
incorporation of new properties such as superconductiv- 
ity or magnetism [Io| . 

The structure of the XYZ Heusler compounds can be 
simply viewed as "stuffed" YZ-zinc-blcndc. Depending 
on the stuffing element X Hcuslers are semiconducting 
or semi- metallic [14|. Materials like ScPtBi are topologi- 



cally similar to HgTe: the inversion of the conduction and 
valence bands occurs due to small electronegativity dif- 
ferences. Since HgTe and ScPtBi are both 2D topological 
insulators, the QSH is also expected in the corresponding 
quantum wells, as e.g. ScPtSb/ScPtBi, in full analogy 
to CdTc/HgTe. We emphasize that the check of parity 
at the time reversal points, as a sufficient condition for 
their topological character, is not possible here because 
of the absence of inversion symmetry Q . 



Heusler compounds are similar to a stuffed diamond, 
correspondingly, it should be possible to find the "high 
Z" equivalent of graphene in a graphite-like structure 
with 18 valence electrons and with inverted bands. In 
this structure type with a lower symmetry compared to 
diamond three-dimensional topological behavior is re- 
alizable. Indeed, the honeycombs KZnP and KHgSb, 
crystallizing in so-called AIB2, M2I11 or ZrBeSi struc- 
ture types, exhibit band ordering similar to the cubic 
CdTe and HgTe. However, in contrast to graphene, these 
compounds have a stron g S OC which leads to a finite 
band gap at the T point [lfil [l6j . The determination of 
the wave function parity at the time-reversal points for 
these materials @ is now possible, since they are centro- 
symmetric. 

In detail, KHgSb, for instance, can be presented as a 
stuffed graphene in the following way: the electropositve 
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K + is stuffed in a honeycomb lattice of [HgSb] ~ |17j . The 
ZrBeSi structure with additional stuffing electropositive 
main group elements is illustrated in Figure [TJ One sees 
that the main group and the transition elements alternate 
within the same layer and between the layers, therefore 
the primitive unit cell consists of two formula units. 

The number of ternary compounds crystallizing in this 
structure is nearly as large as in the Heusler family, how- 
ever due to a weaker mean hybridization (sp 2 versus sp 3 ) 
many of the 18 valence electron compounds exhibit sub- 
stantially smaller gaps or are simply metallic. As we 
show in the following, a certain number of semiconduct- 
ing materials of this family with ordinary and inverted 
band structures can be found. In contrast to HgTe and 
ScPtBi, such materials with inverted structure exhibit 
real band gaps due to their higher structural anisotropy. 
For this reason the honeycomb compounds are three di- 
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FIG. 1. (color online) Structure of KZnP, KHgSb, KCuSe and 
LiAuTe (ZrBeSi structure type) (a) shown in the a-b plane 
and (b) along the c axis. Thus, the ternary semiconductor 
(with 18 valence electrons) can be viewed as stuffed relative 
of graphite with the honeycomb lattice consisting of alternat- 
ing late transition metals (Ag, Au, Zn, Hg) and main group 
elements (Se, Te, P, As, Sb). 



mensional (3D), rather than 2D topological insulators. 
The 3D topological insulators exhibit an insulating en- 
ergy gap in the bulk and gapless states on the surface 
protected by time-reversal symmetry Q. These surface 
states exist only if there is an odd number of massless 
Dirac cones, with a single cone in the simplest case. The 
oddness is provided by the Zi topological invariant of the 
bulk (TEl [la ] , thus any time- reversal invariant perturba- 
tion cannot open an insulating gap at the Dirac point on 
the surface. 

The bulk band gap is defined by interplay of the lattice 
constant, spin-orbit coupling and the difference in elec- 
tronegativities of the honeycomb sublatticc constituents. 
In analogy to a ternary cubic semiconductors, new topo- 
logical insulators can be found among the heavy 8/18 
valence electron relatives of graphene, i. e. in a graphite 
XYZ structure type. Suitable atomic combinations are: 
(i) X = Li, Na, K, Rb, Cs; Y = Zn, Cd, Hg and Z = P, As, 
Sb, Bi, or (ii) X = K, Rb, Cs; Y = Ag, Au and Z = Sc, Te, 
or (iii) X = rare earth, Y = Ni, Pd, Pt and Z = P, As, Sb, 
Bi. Some of these combinations are not yet synthesized, 
some crystallize in different forms such as variants of the 
Ci^Sb structure type (they will be the subject of another 
publication) and some, especially the rare earth contain- 
ing compounds, arc metallic. The known examples are 
KZnAs, KZnSb, KZnP, KHgAs, KHgSb, RbZnP, RbZ- 
nAs, RbZnSb, NaAuTe, KCuSe, KCuTe, KAuTe, and 
RbAuTe. Compounds such as LiAuSc, LiAuTe, CsAuTc, 
KHgBi, and CsHgBi are likely to be synthesized and we 
consider them in the present study as well. As mentioned 
above, the electronic structure of these materials is simi- 
lar to their cubic analogues. The alkaline ions Li + , Na + , 
K+, Rb+, and Cs+ "stuff" the graphite type YZ pla- 
nar sublattice. Since these 18-clcctron compounds form 
such closed-shell structures, they are all non-magnetic 
and semiconducting. The difference is that in the case of 
binary semiconductors or Cu, Heuslers, the bonds within 
the YZ tetrahedrons are of sp 3 or sd 3 type, whereas 
in planar graphite-type layers the cr-type bonding oc- 
curs between the sp 2 or sd orbitals. The remaining p- 
orbitals provide the 7r-type bonding interaction, similar 



to graphite. 

The search for the topological character of the pro- 
posed materials is based on ab-initio calculation of the 
electronic structure (for details see the supplemental). 
An important peculiarity of hexagonal systems is em- 
phasized here: their unit cell consists of two formula 
units which results in a doubling of the correspond- 
ing bands. More details are given in the supplemental 
(Fig. 1) where it is seen that for KZnP (KHgSb) at the 
r point the doubled s-bands are split by about lOmcV. 
Thus for the weak coupling of the nearest honeycomb 
planes these doubled bands can become nearly indistin- 
guishable. In the following we show that such doubling 
of the unit cell leads to a change of topology. Fortu- 
nately, the centro-symmetric space group (194) of the 
honeycomb-type compounds allows to make use of the 
parity eigenvalues @, H[. All relevant properties, i.e. 
the wave function parity in time-reversal symmetric k- 
points (r(0, 0, 0), M (tt, 0, 0), L{ir, 0, tt), A(0, 0, tt)), the Z 2 
invariant, the average nuclear charge (Z), the band gap 
width E g and corresponding lattice parameters are listed 
in Table I for the compounds studied here. 



TABLE I. Wave mnctionparities at the time-reversal points, 
Z-2 topological invariant [3], average nuclear charge (Z), band 
gap E s and optimized lattice parameters. 
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From previous studies on binary and ternary cubic 



FIG. 2. Band structure of LiAuSe (a, b) and KHgSb (c, d) 
calculated without (a, c) and with (b, d) SOC. 



semiconductors [lfj it follows that it is more probable 
to find the Zi = 1 topological insulator among heav- 
ier compounds (with stronger SOC). Indeed, their bands 
splitting scales roughly with the average nuclear charge 
(Z) = l/NJ2f=i Zi where N = 2 for binaries and TV = 3 
for ternaries. This parameter sorts cubic systems almost 
along a straight line However Table I clearly illus- 
trates that this does not hold for the semiconductors of 
the ZrBeSi structure type. Indeed, only Lithium com- 
pounds within the Zn and Hg group show the expected 
trend. The compounds with Zn are topologically trivial 
whereas those with Hg are non-trivial insulators. Com- 
pounds with heavier alkaline metals (Na, K) are all triv- 
ial independently of whether or not they contain Zn or 
Hg. Among the Ag- and Au-containing compounds more 
non-trivial systems are found, however the correlation 
between band inversion and (Z), as in the cubic semi- 
conductors, is absent. For example, the topological in- 
sulator LiAgSc corresponds to (Z) = 28, whereas for the 
topologically trivial KHgSb system (Z) = 50. 

In Figure[2]the band structure of LiAuSe (upper panel) 
is compared with KHgSb (lower panel) calculated at 
time-reversal symmetric points with (right) and without 
(left) SOC. It follows that both compounds are semi- 
metals with degeneracies at the Y and A symmetry points 
if SOC is omitted. The degenerate p x and p y states me- 
diate a- type Sb/Se and Hg/Au bonding. The lower-lying 
Hg/Au bands are of s type, similar to HgTe. Inclusion of 
SOC opens a band gap at the Fermi energy at the F and 
A points leading to the typical dips in the band structure 
for both compounds [llj]. The resulting parities for both 
materials listed in Table I match only at the M point, 
whereas at T, A, and L they differ, leading to a trivial 
state in KHgSb and topological insulator in LiAuSe. 

To understand the mechanisms of the band inversion 
and parity change in details one can track the band 
structure evolution near the T point by starting from 
the simple atomic energy levels and subsequently in- 
troducing chemical bonding (step I), crystal field (II), 
and SOC (III). The resulting changes are schematically 
shown in Fig. [3] As an example we took the KZnP sys- 



FIG. 3. (Color online) The structure of the atomic orbitals. 
(a) Evolution of atomic orbitals at the Y point in KZnP sys- 
tem, which has no inverted bands formed by Zn s and P 
Px.y.z-orbitals. Steps (I-III) represent the effect of turning on 
chemical bonding (I), crystal-field splitting (II) and SOC (III). 
The blue dashed line marks the Fermi energy. Similar scheme 
for LiAuSe (b) and KHgSb (c), (d), (e) and (f) represent the 
schematic band structure for KZnP, LiAuSe and KHgSb near 
the r point, respectively. The black curves mark the orbitals 
with mainly Zn (or Au and Hg) s-character, red and green 
curves are the rest of light and heavy hole-like bands. 

tern. The low-lying 3s-orbital of P(3s 2 3p 3 ) can be ne- 
glected and the consideration can be restricted to the 
s-orbitals of K(4s) and Zn(3d 10 4s 2 ) and the p x ,y,z or- 
bitals of P(3s 2 3p 3 ). In step I the chemical bonding is 
introduced. It is convenient to make use of the inver- 
sion symmetry and recombine the orbitals according to 
their parities. K and Zn give two even and two odd s- 
orbitals |-Sk,±)) \Szn,±)', P gives three odd and three 
even p-orbitals |Pp )Xyz ,±), where "±" are the parity la- 
bels. In step II the crystal field is switched on which splits 
the |Pp, X yz,±) into |Pp, xy ,±) and |Pp, z ,±) according to 
the hexagonal symmetry. In step III the SOC is turned 
on. This leads to a splitting between |Pp )X +iy,±,4.(t)) and 
Pp,z,±,t(4-)) orbitals. The Fermi level falls into the mid- 
dle of |Pp,x+iy,-,i) (or |Pp, x -iy,-,t)) and |Pp, x+iyj+lt ) 
(or IPp^-iy^j,)) as shown in Figure [3](a) . No inversion 
occurs between the s-orbital of K (or Zn) and p-orbital 
of P, similar to the CdTe case. However, this situation 
changes for the LiAuSe and KHgSb compounds. Due 
to the very delocalized character of Au (Hg) d-orbitals 
the s-orbital of Li (K) is pulled down below the p-orbital 
of P, which leads to band inversion similar to the HgTe 
case [lfj. The corresponding evolution for LiAuSe and 
KHgSb is shown in Figures [3](b) and (c) . In addition, 
|Psb, x +iy,-,4.) (or |Psb, x -iy,-,t)) and |Psb,x+i y ,+,t) ( or 
|Psb, x -iy,+,j,)) are inverted because of the strong SOC 
in KHgSb. We emphasize that since there are two band 
inversions occurring in KHgSb it becomes topologically 
trivial. In contrast, LiAuSe exhibits inversion only be- 
tween the s- and p-orbitals since its SOC is too weak to 
invert |Psb, x +i y ,-,i) (or |Psb, x -i y ,-,t)) and |Psb,x+iy,+,t) 
(or |Psb, x -iy,+,4.))- Thus we conclude that LiAuSe is 
topologically non-trivial. The corresponding band struc- 
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(a)KZnP 


(b) KHgSb / 








FIG. 4. (color online) The band structure of the semi-infinite 
(vacuum/solid-interface) corresponding to (a) trivial KZnP, 
(b) non-trivial LiAuSe and (c) trivial KHgSb insulators, cal- 
culated along the K-T-M directions of the Brillouin zone. In 
contrast to KZnP and KHgSb, LiAuSe exhibits the surface 
states seen in the bulk energy gap (zoomed in the inset (d)). 



tures around the V point for KZnP, LiAuSe and KHgSb 
are shown schematically in Figurcs[3](d-f). It is seen that 
in KZnP no band inversion occurs, in LiAuSe it occurs 
once and KHgSb twice. 

For the final test of the topological/trivial character we 
directly calculate the surface states of the proposed sys- 
tems. The electronic structure of the surface can be ade- 
quately described by the vacuum/solid interface within a 
unit cell which is sufficiently large along a certain direc- 
tion. By utilizing the so-called decimation technique [2l[ 
such a unit cell can be extended to infinity which ensures 



correct boundary conditions. Such calculations typically 
require a formalism using fast decaying basis functions. 
In the present work we apply the Wannier functions ap- 
proach [22| • As it follows from the band structures shown 
in Figure 0] only LiAuSe exhibits gapless surface states 
seen as a single Dirac cone within the bulk energy gap 
at the r point. The estimate of its Fermi velocity gives 
about 1.8 x 10 5 m/s which is smaller than that for Bi2Sc3. 
In contrast, KZnP and KHgSb remain insulators at the 
surface. Thus the surface state calculation agrees with 
the bulk parity analysis and conclusively confirms the 
existence of the topologically non-trivial materials within 
proposed honeycomb type structure. 

In conclusion we emphasize that topologically non- 
trivial systems can be found in the proposed class of 
honeycomb structure semiconductors. The interplay of 
mechanisms which are responsible for the topologically 
trivial or non-trivial character in these systems differs 
from the cubic semiconductors studied earlier. In partic- 
ular it is shown that the strong SOC and weak inter-layer 
coupling causes a double inversion which in turn makes 
the compound trivial. In contrast to the topologically 
non-trivial cubic systems which exhibit a zero band gap 
in the bulk, the topologically non-trivial hexagonal mate- 
rials provide the "natural" 3D topological materials with 
a real bulk band gap and gapless states at the surface. 
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Calculation of the electronic structure is performed by 
the plane-wave based pseudopotential BSTATE (Beijing- 
Simulation Tool of Atomic TEchnology) package [1| . The 
exchange-correlation potential is treated within the Gen- 
eralized Gradient Approximation (GGA) Q. The lattice 
parameters listed in present work are obtained by doing 
the ab-initio geometry optimization. 

The topological character is calculated by following 
Ref. Q. The central quantity here is the parity prod- 
uct Si = Yl m=1 £2™ (Pi), where N is the number of filled 
bands and ^mij'i) is the parity eigenvalue of the 2m-th 
occupied band at the time reversal point T^. The Z2 in- 
variant, 1/0 = or 1, is then obtained from the product 

(-ir° = n^- 

Figure [1] illustrates the analogy between cubic Cu 
compounds and hexagonal semiconductors, by compar- 



ing the band structures of CdTe and HgTe with hexag- 
onal ternaries KZnP and KHgSb. Clear fingerprints can 
immediately be identified: both CdTe and KZnP exhibit 
a direct gap at the T point formed by splitting of the 
conduction (s-type) and the valence (p-type) bands. By 
going from CdTe to HgTe and from KZnP to KHgP the 
bands at the T point invert in a similar manner: the s- 
type shifts below the p-type band. However, in contrast 
to HgTe and other Cib Heuslers, the hexagonal symme- 
try with varying c/a ratios allows to lift the degeneracy 
of the p-typc bands and opens a finite band gap. As we 
mentioned this suggests KHgSb and the related materials 
as candidates for a 3D topological insulators, similar to 
Bi2Se3. It follows from the fiat bands along the K — H di- 
rection that the material is quasi two-dimensional which 
leads to a strong anisotropy of its transport properties. 
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FIG. 1. (color online) Comparison of the band structures of 
cubic binaries (CdTe, HgTe) with hexagonal (KZnP, KHgSb) 
ternary compounds. In both CdTe and KZnP the valence and 
conduction bands are formed mainly by the p (red) and the s 
(blue) states, respectively. 



